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Abstract 

A new mathematical model for the Black-Scholes equation is proposed to forecast 
option prices. This model includes new interval for the price of the underlying stock 
as well as new initial and boundary conditions. Conventional notions of maturity 
time and strike prices are not used. The Black-Scholes equation is solved as a 
parabolic equation with the reversed time, which is an ill-posed problem. Thus, a 
regularization method is used to solve it. This idea is verified on real market data 
for twenty liquid options. A trading strategy is proposed. This strategy indicates 
that our method is profitable on at least those twenty options. We conjecture that 
our method might lead to significant profits of those financial institutions which 
trade large amounts of options. We caution, however, that detailed further studies 
are necessary to verify this conjecture. 

Key Words: Black-Scholes equation, new mathematical model, new initial and 
boundary conditions, testing on real market data, parabolic equation with the reversed 
time, Ill-Posed problem, regularization method 


1 Introduction 

The Black-Scholes equation is solved forwards in time to forecast prices of stock options. 
Since this is an ill-posed problem, a regularization method is used. Uniqueness, stability 
and convergence theorems for this method are formulated. For each individual option, 
we use its past history as an input data. We propose a new mathematical model for the 
interval of prices of the underlying stock, as well as for initial and boundary conditions 
on this interval for the Black-Scholes equation. The conventional notions of strike prices 
and maturity time are not used. Our model does not make distinctions between put and 
call options. Based on our results for market data, we conjecture that our methodology 
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might result in significant profits if trading a large number of options, as it is done in 
some large hedge funds. We caution, however, that additional studies of a substantial 
number of options are necessary to verify this conjecture. 

To verify the validity of our model, we use the real market data for twenty liquid 
traded options. We believe that this is the best possible way for the verification. We 
select options randomly among liquid ones. Those liquid options were taken selected at 
http://finance.yahoo.com/options/lists/, Market prices, implied volatility and prices of 
the underlying stock for our selected options were taken from the Bloomberg terminal 

http://www.bloomberg.com, The single condition we impose when selecting these 
options, is that these options should be daily traded. Based on our technique, we propose 
here a certain trading strategy. This strategy shows that we are profitable on seventeen 
(17) out of twenty (20) options. 

Our mathematical model allows one to forecast option prices of daily traded options 
from the current time event to two next ones, i.e. for a short time period. As to the 
time distance between two events, it depends on time units which one is using. In our 
calculations one time unit is one trading day. Thus, we forecast prices for “tomorrow” and 
the “day after tomorrow” having knowledge of past prices of “the day before yesterday”, 
“yesterday” and “today”. For our convenience, we use only last prices in each particular 
day, i.e. prices which were in place just before the closure of the market. But one unit 
can be one minute, one hour, one week, etc., and also last prices can be replaced with 
prices at other moments of time. An optimization of the latter is a technical issue which 
can be addressed later. 

In the conventional approach to the Black-Scholes equation, one solves this equation 
backwards in time t G (0, T ), having an initial condition at the maturity time t = T. 
The initial condition at t = T includes the strike price K, see Hull (2000) and Wilmott, 
Howison and Dewyne (1997). However, there is a major drawback in this approach. In¬ 
deed, T is usually a few months. But it is obviously impossible to forecast the value of 
the volatility for such a long time period with a reasonable accuracy. On the other hand, 
the solution of the Black-Scholes equation critically depends on the volatility coefficient. 

The above led us to believe that it is more natural to use the Black-Scholes equation to 
forecast prices of stock options for short time periods. However, the first major obstacle 
in this direction is that the problem of solving the Black-Scholes equation forwards in 
time is ill-posed. Namely, its solution usually does not exists and even if it exists, then 
it is unstable with respect to the input data. Therefore, a regularization method should 
be used. The second major obstacle is that some crucial input parameters for solving 
that equation forwards in time are unknown. These parameters are: the interval of prices 
of the underlying stock, the initial condition, two boundary conditions and the volatility 
coefficient. 

We positively address here the following two questions: 

Question 1. Is it possible to forecast the option price for a rather short time period 
using the Black-Scholes equation and market data? 

Question 2. If this is possible, then can one be profitable using this forecast and a 
trading strategy? 

There are four main questions which need to be addressed to answer the first question: 

1. What is the interval for the prices of the underlying stock? 

2. What are boundary and initial conditions on this interval? 
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3. What are the values of the volatility coefficient in the future? 

4. How to solve the Black-Scholes equation forwards in time? 

We address questions 1-3 in our new mathematical model. To address the fourth 
question, we use a regularization algorithm which was developed in Klibanov (2014). We 
formulate theorems about stability and convergence of this method. They were proven in 
Klibanov (2014) for a general parabolic equation of the second order. The key method of 
proofs in these references is the method of Carleman estimates. 

We now briefly explain why the solution of the Black-Scholes equation in the forward 
direction of time is an ill-posed problem. Let r = const. > 0 and the function / G 
L 2 (0,7r). Here is a well known example demonstrating the ill-posedness of problems for 
parabolic equations with the reversed time. Consider the following problem for the heat 
equation with the reversed time, 

u t + u xx = 0, (x, t) e (0,7r) x (0, r), 

u (x, 0) = / (x ), 
u (0, t) = u ( 7 r, t ) = 0. 

It is well known that this problem has no more than one solution, see, e.g. Klibanov 
(2014) and Lavrentiev, Romanov and Shishatskii (1986). The unique solution of this 
problem, if it exists, is 

OO 

u (x, t) = f n sin (nx) e nH , 

71=1 

where {f n }™ = i are Fourier coefficients of the function f (x) with respect to sin (nx). 
Consider the L 2 (0, 7 r) norm of the function u (x, t ), 

OO 

MM)iiL ( o ,„) = £* 2 ” 2 ‘- (i-i) 

71=1 

Hence, if the solution of this problem exists, then squares of Fourier coefficients must 
decay exponentially with respect to n. Hence, the solution of this problem exists for only 
a narrow set of functions /. Also, it is clear from (11.11) that even if this series is truncated, 
still small fluctuations of can lead large variations of the function u. This manifests 
the instability of the above problem. Furthermore, it follows from (11.11) that the larger r 
is, the more unstable this problem is. Hence, in order to obtain a rather good accuracy, 
any regularization method should work only on a short time interval. Thus, a rather 
accurate forecast of option prices via the Black-Scholes equation might occur only for a 
rather short time period. The latter is exactly what we are doing here. 

In section 2 we present our mathematical model. In section 3 we describe our numerical 
method. In section 4 we show our results for market data. We summarize our results in 
section 5. 
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2 The new mathematical model 

Let s be the stock price, t be time, a (t) be the volatility of the option. In onr particular 
case we use only the implied volatility listed on the market data of http://www.bloomberg.com. 
However, more sophisticated models for the volatility can also be used in our model. Let 
r > 0 be the unit of time for which we want to forecast the option price. In our partic¬ 
ular case r is one trading day since we forecast option prices for “tomorrow” having the 
knowledge of these prices, as well as of other parameters for “today”, ’’yesterday” and 
“the day before yesterday”. Since there are 255 trading days per year, then in our case 
r = 1/255. However, our model can work with the case when r is any unit of time: one 
hour, one minute, etc. In any case we assume below that r G (0,1/4). Let ’’today” be 
t — 0, “tomorrow” be t — r, the day after “tomorrow” be t = 2 r, “yesterday” be t = — r 
and “the day before yesterday” be t — —2 r. 

We forecast only the last price of the option, i.e. the price on which one item of 
that option was either bought or sold in the last trading event of a trading day for this 
particular option. Let Ub (t) and u a (t) be respectively the bid and ask prices of the option 
at the moment of time t. Let Sb (t) and s a (t) be respectively the bid and ask prices of the 
stock at the moment of time t. It is well known that Ub (t) < u a (t) and that Sb (t) < s a (t ). 
Thus, the market data which we need for our model are listed in Table 1. These data are 
available at http://www.bloomberg.com, 

Table 1. The data we need to forecast the last option price for the next two 

days t — T and t = 2 r. 







Ub it) 

U a (t) 
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t = — 2r, —r, 0 

t = — 2t, —r, 0 

t = 0 

t = 0 


First, having discrete values of functions Ub (t ), u a (£), a (t) for three moments of time 
listed in this table, we interpolate these functions between these three points using the 
standard quadratic interpolation. Thus, we obtain approximate values of these functions 
for t G (—2t, 0) as quadratic polynomials. Next, we extrapolate functions Ub(t),u a (t) 
for t G (0, 2r) as those quadratic polynomials. Since r is small, then it is reasonable to 
assume that functions Ub (t) , u a (t) , a (t) are approximated rather well on the time interval 
t G (—2r, 2r) . Thus, we obtain functions 

u b (t) , u a (t) , a (t) for t G (0, 2r). 

Let u = u(s, t ) be the price of one item of the stock option. Denote Sb = Sb (0) , s a = 
s a (0). Then Sb < s a - Consider the numbers Ub = u (s&, 0) and u a = u (, s a , 0). Let 


f(s) 


U b ~ Ug 

Sa 


li a Sb 'U'bSa 

• s H- 

Sb S a 


( 2 . 1 ) 


be the linear interpolation between Ub and u a on the interval s G (sb, s a ). Hence, / (s&) = 

Ubi f (®a) U a . 

The simplest form of the Black-Scholes equation is 


Lu : = 


u t + 


cr 


(t) 


S 2 u s 


2 


0, 


( 2 . 2 ) 
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where L is the partial differential operator of the Black-Scholes equation, see Hull (2000) 
and Wilmott, Howison and Dewyne (1997). We solve this equation on the interval of 
stock prices s G (s b , s a ) and for times t G (0, 2r). Thus, we impose the following initial 
condition at t — 0 and boundary conditions at s = s b , s a : 


u (s, 0) = /(s),sG (s b ,s a ), 
u (s b , t ) = u b (t ), u (s a , t ) = (t) , for t E (0, 2r). 


(2.3) 

(2.4) 


Denote Q 2t = {(s,i):sG (s b ,s a ),t G (0, 2r)}. In this paper, we computationally 
solve the following problem: 

Problem. For (s, t ) G Q 2t find the solution u (s, t) of equation \2.2\) satisfying the 
initial condition 112. ,9j) and boundary conditions 

This problem as well as the above interpolations and extrapolations form our new 
mathematical model. Theorem 1 claims uniqueness of the solution of this problem. 
This theorem follows immediately from Klibanov (2014) and Lavrentiev, Romanov and 
Shishatskii (1986). Below H 2,1 (Q 2t ) an d H 2 (Q 2t ) are standard Sobolev spaces of real 
valued functions. 

Theorem 1. The problem \2.2\)-\2.f\) has at most one solution u G H 2,1 (Q 2t ) . 


3 Numerical method for the problem (I2.2l )- ([2T4T ) 


As it was pointed out in Introduction, the problem fj2.2h - fj2.4jl is ill-posed, since we are 
trying to solve the Black-Scholes equation forwards rather than backwards in time,. The 
ill-posedness means here that the existence of the solution is not guaranteed. In addition, 
the solution, even if it exists, is unstable with respect to small fluctuations of initial and 
boundary conditions, see, e.g. the book of Tikhonov, Goncharsky, Stepanov, and Yagola 
(1995) for the theory of Ill-Posed Problems. Therefore, we use the regularization method 
of section 5 of Klibanov (2014). In simple terms, we find such an approximate solution 
of this Problem, which satisfies conditions fj2.2h - fj2.4jl in the best way in the least squares 
sense. 

Consider the following function F (s, t) 


p ^ ^_ n b (fi) u a (f) s p n a (t) s b u b (i) s a 


$b Sa Sjj S a 

Then F G H 2 (Q 2t ) . It follows from (12 .1 jl . (12.311 . (12.4h and (13.ip that 
F (s, 0) = / (s), F (s b , t) = u b (' t ), F (s a , t ) = u a (t). 


(3.1) 


(3.2) 


Following section 5 of Klibanov (2014), consider the following Tikhonov-like functional 


J a (u) = j ( Lu ) 2 dxdt + a ||u — F\\h 2 (q 2t ) , (3.3) 

Q2t 

where a G (0,1) is the regularization parameter. Note that in the conventional case of 
linear ill-posed problems Tikhonov functional is generated by a bounded linear operator, 
see, e.g. Ivanov, Vasin and Tanana (2002). However, in our case L is an unbounded 
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differential operator L : H 2,1 (Q 2t ) —>■ L 2 ( Q 2t ), in which case H 2,1 (Q 2t ) is considered as 
a dense linear set in the space L 2 (Q 2 t) ■ We consider the following minimization problem: 

Minimization Problem. Minimize the functional J a (u) in \S.S 1) . subject to the 
initial and boundary conditions 112.Si) . \2.j\ ). 

To solve this minimization problem computationally, we have written partial deriva¬ 
tives in (13.31) via finite differences. In particular, we have obtained a finite difference grid 
covering the rectangle Q 2r . Next, we have minimized J a (u) with respect to the values of 
the function u(s,t ) at grid points, using the conjugate gradient method. The starting 
point of this method was u = 0. 

The existence and uniqueness of the minimizer of functional (13.31) supplied by boundary 
conditions (12. 3[) , (12.4p is ensured by Theorem 2. This theorem follows immediately from 
Theorem 5.3 of Klibanov (2014), where the general parabolic equation of the second order 
was considered. 

Theorem 2. In \3.3 1) . let F be the function defined in \3.2 1) . Then for any a G (0,1) 
there exists a unique minimizer u a G H 2 (Q 2t ) of the functional J a (u) satisfying condi¬ 
tions $2.3\), Furthermore, there exists a constant C = C (Q 2t , a) > 0 depending 

only on listed parameters such that 

C 

\\ u »\\h^(Q 2 T ) - fjff W F \\h^q 2t ) ■ 

To formulate the convergence result for minimizers u a , we need to assume that there 
exists the “ideal” exact solution of our ill-posed problem, i.e. solution with the ideal 
noiseless data. Such an assumption is one of most important building blocks of the 
Tikhonov regularization theory, see, e.g. Tikhonov, Goncharsky, Stepanov, and Yagola 
(1995). Thus, we assume that there exists the exact solution u* (s,t) G H 2 (Q- 2t ) of 
equation (12.2j) with exact initial condition f* (s) G H 2 ( Sb , s a ) in (12.31) and exact boundary 
conditions ul (t) i u * a ( t ) G H 2 (0,2r). We also assume that there exists a function F* G 
H 2 ( Q 2t ) such that 

F* (s, 0) = f* (s ), F* (s b , t ) = u* b (; t ), F* (s a , t ) = < (t ), (3.4) 

which is similar with (13.21) . Let 6 G (0,1) be a sufficiently small number. We assume that 

ll-^ 1 - -^*H^ 2 (Q2r) - ( 3 - 5 ) 

Hence, it follows from (13.21) . (13.41) and (13.51) that the number 6 characterizes the level of the 
error in our data / (s), u b (t) , u a (t) as compared with the exact data f* (s), ul (t) , u* (t). 

In Theorem 3 we estimate the convergence rate of our minimizers u a to the exact 
solution u*, assuming that 8 —> 0. Theorem 3 follows immediately from Theorem 5.4 
of Klibanov (2014), which was proven for a general parabolic operator of the second 
order. As it is always done in the regularization theory, we choose in Theorem 3 the 
regularization parameter a = a (8) depending on the level of the error in the data. Note 
that assumption (13.51) of the small perturbation in the data seems to be close to the 
reality. Indeed, although we do not assume that the function f* (s) is linear and functions 
ul ( t ), w* ( t ) are quadratic ones, as our functions / (s), Ub ( t ), u a ( t ) are, still since intervals 
(sb, s a ) and (0, 2r) are small, then it is reasonable to assume, as we do, that the function 
/ (s) is a linear and functions Ub ( t ), u a ( t ) are quadratic. 
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Theorem 3. Assume that conditions ([ 3.f\) and 113.5 ]) hold. Let a G (0,1) be a number, 
the function a ( t ) G C 1 [0, a] and there exist constants a 0 , cr x > 0, cr 0 < ay srtch that a (t) G 
[ao,cri] for t G [0, a] . Choose the regularization parameter a = a (S) = S 2/3 , where (3 = 


const. G (0,1). Then there exists a sufficiently small number r o = ro ( (x 0 , a i, HcrU^i^o a 


(0, a] depending only on listed parameters such that if r 
convergence rates hold 


G (0, r 0 ), then the following 


II &“«(*) - d * u *\\ l 2 {q t ) + ll M «(5) - u *\\l 2{Qt) < b {} + IMI h*(q 2t )) 57 » ( 3 - 6 ) 

K(*) ( s ) 2r ) — ( s -> 2r ) IIl 2 (s6,s q ) ^ { 1 + W u *\\hhq 2t )) ’ ( 3 - 7 ) 

V ln ) 

where the number 7 = ((3 ln2) /In 4 and the constant B = B (a,Q 2 T ,cr) > 0 depends only 
on listed parameters. 

It is clear from this theorem that the accuracy estimate (13.61) in a smaller domain 
Qt C Q 2t is of the Holder type. On the other hand the estimate (13.71) on the upper side 
of the rectangle Q 2t is of the logarithmic type. Clearly, the accuracy guaranteed by (13.61) 
is better than the accuracy guaranteed by (13.71) . Hence, one can expect to obtain more 
accurate computational results in Q r , as compared with those in Q 2r . This is exactly the 
reason why we have chosen to solve the problem (I2.2l) - (l2.4p in the larger domain Q 2t . It is 
also clear from Theorem 3 why we have chosen the number r > 0 to be sufficiently small. 


4 Results for the market data 

In our computations we have always used r = 1/255, as mentioned in section 2. In other 
words, we forecasted option prices for two days ahead, i.e. for t G (0,2r). To make 
sure that we obtain accurate results at least for the computationally simulated data, we 
have solved first the above Minimization Problem for synthetic data. To simulate these 
data computationally, we have solved equation (12.21) downwards in time with boundary 
conditions (12.41) and the initial condition at t = 2r. This way we have obtained the 
function u (s, 0) = f S i m ( s ). Next we have solved equation (12.21) with boundary conditions 
(12.4p and the initial condition f sim (s) upwards in time by the above described method. 
Next, we have compared the resulting solution with the computationally simulated one. 
Results were quite accurate ones. I 11 this study of computationally simulated data we 
have found that a = 0.01 is the optimal value of the regularization parameter and we 
have used a = 0.01 in all follow up computations. 

Next, we have used the market data which we took at the Bloomberg terminal 

http://www.bloomberg.com, We have compared our results with the true last prices 
ui (r) and ui (2r). Let s m = (s?, + s a ) /2 be the mid point between bid and ask “today’s” 
prices of the underlying stock. First, we have calculated the minimizer u a (s,t) via the 
above minimization procedure, where ( s,t ) G Q 2t . This was done for every option for all 
days of its existence. 

Those cases when u a ( s m , r) was between our extrapolated values of bid Ub (r) and 
ask u a (r ) prices, i.e. within the bid/ask spread, were not of any interest to us and we 
have not traded options in those days. However, we were only interested in those cases 
when our predicted price was larger than our extrapolated ask price for at least $0.02, i.e. 















u a (s m , r) > u a (r) + 0.03. Indeed, bid and ask prices are close to each other, which makes 
the case uj, (r) < u a (s m , r) < u a (r) not interesting for trading. On the other hand, the 
cut-off value of $0.02 was chosen due to our computational experience. The same was for 
u a ( s m , 2r), see our trading strategy below. In all cases of trading strategy listed below 
we buy and sell options just before the market closure. 

Trading Strategy: 

1. Suppose that u a ( s m , r) > u a (r) + $0.03 and also that u a ( s m , 2r) > u a (2r) + $0.03. 
Then we buy two (2) items of that stock option at t = 0, i.e. “today”. Next, we sell 
one item at t — r and sell the second item at t — 2r.In doing so, we do not make 
any forecast at t — r. Next, we apply our forecasting procedure at the day t = 2r 
as above and repeat the use of our trading strategy. 

2. Suppose that u a (s m , r) > u a (r) + $0.03 but u a (■ s m , 2r) < u a (2r) + $0.03. Then we 
buy one (1) item of that stock option at t — 0. Next, we sell that item at t — r. 
Next, we apply our forecasting procedure at the day t = t as above and repeat the 
use of our trading strategy. 

3. Suppose that u a (s m , r) < u a (r) + $0.03, but u a (s m , 2r) > u a (2r) + $0.03. Then 
we buy one (1) item of that stock option at t — 0. Next, we sell that item at t = 2r. 
Next, we apply our forecasting procedure at the day t = 2r as above and repeat the 
use of our trading strategy. However, we do not apply our forecasting procedure at 
the day t — r. 

4. Suppose that u a (s m , r) < u a (r) + $0.03 and also that u a ( s m , 2r) < u a (2r) + $0.03. 
Then we neither buy nor sell this option “today”. Note that a particular case of 
this is the situation when both forecasted prices for t — r and t = 2r are within 
the bid/ask spread, i.e. when (r) < u a (s m ,r ) < u a (r ) and also Ub (2r) < 
u Q (. s m , 2r) < u a (2r). 

We have evaluated twenty (20) liquid options, which were selected randomly among 
those options which are daily traded (see Introduction). Their short codes and our num¬ 
bers for them are given in Table 2. 

Table 2. Short option codes. “C” and “P” mean call and put options 
respectively. “Number of days” means the total number of days an option 

was evaluated by our procedure. 
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Option Number 

Option Short Code 

Number of days 

1 

INTO US 01/17/15 C25 Equity 

181 

2 

WFC US 11/22/14 C50 Equity 

52 

3 

PFE US 11/22/14 C29 Equity 

68 

4 

YHOO US 01/17/15 C50 Equity 

267 

5 

AIG US 11/22/14 C55 Equity 

155 

6 

AAPL LIS 01/17/15 C80 Equity 

130 

7 

SBUX US 12/20/14 C80 Equity 

30 

8 

AAL US 01/17/15 C40 Equity 

195 

9 

QCOM US 12/20/14 C80 Equity 

31 

10 

MSFT US 01/17/15 C45 Equity 

267 

11 

QQQ LIS 11/22/14 C106 Equity 

45 

12 

MRK US 11/22/14 C60 Equity 

58 

13 

DDD US 01/17/15 C45 Equity 

119 

14 

WMB US 02/20/15 C48 Equity 

43 

15 

MSFT US 03/20/15 C50 Equity 

48 

16 

SPY US 02/20/15 P205 Equity 

141 

17 

1WM LIS 03/20/15 PI 10 Equity 

115 

18 

EEM US 02/20/15 P39 Equity 

108 

19 

YHOO US 03/20/15 C55 Equity 

76 

20 

HYG US 03/20/15 P86 Equity 

48 


Table 3 shows total profit/loss for each option resulting from our price forecast and 
the follow up application of our trading strategy. Losses are with the ” sign. We did 
not invest any real money. Rather, for each option, we pretended that we buy and then 
sell the next day only one item and only at those days which were recommended by our 
trading strategy. To get results of this table, we have compared u a (s m ,r) with the true 
last price of the next day ui (r). 


Table 3. Total gain/loss for each option of Table 2 during the evaluation 
period. The last row shows the total profit for all options. 
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Option Number 

Profit/loss 

1 

2.69 

2 

0.51 

3 

0.51 

4 

2.93 

5 

1.79 

6 

4.72 

7 

0.13 

8 

2.65 

9 

0.65 

10 

1.74 

11 

0.008 

12 

0.9 

13 

-1.66 

14 

-0.63 

15 

0.37 

16 

0.98 

17 

-2.26 

18 

2.9 

19 

1.06 

20 

0.77 

Total 

20.8 


The next interesting question is: How accurate are we in our forecast? Table 4 gives 
the average relative error of our forecast for each option. For each option, we have 
calculated the relative error for those days when we were selling that option by the above 
strategy. That relative error was 


\u a ( s m ,r ) -ui(t) | 
ui (r) 

Next, we took the average value of numbers (14.1|) for each option of Table 2. 


(4.1) 


Table 4. The average relative accuracy of our forecast. The last row shows 
the total average relative accuracy over all options. 
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Option Number 

Relative Error 

1 

0.179935747 

2 

0.282359905 

3 

0.317038779 

4 

0.197047078 

5 

0.201012726 

6 

0.217644086 

7 

0.295314034 

8 

0.243892569 

9 

0.270193217 

10 

0.237299767 

11 

0.243110472 

12 

0.172772735 

13 

0.173193171 

14 

0.371006086 

15 

0.135316766 

16 

0.214612849 

17 

0.27188249 

18 

0.258101445 

19 

0.22011074 

20 

0.245692644 

Average for All 

0.237376865 


Table 4 shows that the average relative error (14.11) is rather large: 23.7%. We have 
observed that we often overestimated option prices. Nevertheless our trading strategy 
enables us to be profitable at least on twenty options of Table 2. 

5 Summary and conclusions 

We have proposed a new mathematical model for the Black-Scholes equation. Instead of 
traditionally solving this equation backwards in time, which is a well posed problem, we 
solve it forwards in time, which is an ill-posed problem. The theory, which was previously 
developed in Klibanov (2014), guarantees the existence and uniqueness of regularized 
solution of this problem as the convergence of regularized solutions to the correct one as 
the level of the noise in the data 5 and the regularization parameter a (5) tend to zero 
(Theorems 2,3). 

The best way to verify a mathematical model is to apply it to the real market data. 
And this is what we have done here. We have chosen randomly twenty (20) liquid options. 
Next, we have presented forecasts of option prices for next two days for them. We have also 
developed a simple trading strategy, which is using our forecast to buy and sell options. 
We pretended that we buy and sell those options (the real money were not invested). The 
main conclusion, which can be drawn, is that even though our technique almost always 
overestimates option prices, still we got profits in seventeen (17) of above options and we 
got losses in three (3) of them. Furthermore, the total result is that, summing up these 
profits and losses, we are still profitable on our randomly selected twenty options. 
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We point out that we did not take into account the transaction cost. We conjecture 
that accounting for this cost would increase the above cut-off value of $0.02, while leaving 
the conclusion about the profit intact. 

Even though our profits are small, this should not be discouraging. Indeed, we pre¬ 
tended to buy and then sell the next day only one-two items of an option at a time, and 
this was done only for twenty (20) options. However, large financial institutions buy and 
sell a large number of options daily. Hence, we conjecture that such large trading might 
indeed be profitable if using the technique of this paper. Still, we caution again that 
further studies of a much larger number of options should be performed to figure out how 
correct our conjecture is. However, we do not have enough resources to conduct such 
studies. 

We believe that it is possible to refine our results and probably to increase potential 
profits. We are naturally concerned about the way of refining the accuracy of our predicted 
prices, see Table 4. One of such ways is to obtain more accurate values of the volatility 
coefficient, as compared with the implied volatility we have used here. Next, one should 
solve the problem fl2.2p - fl2.4p with this more accurate volatility coefficient. Most likely, 
the volatility depends on both stock price s and time t, a = a (s, t). However, to compute 
the volatility, one needs to solve a very difficult coefficient inverse problem for the Black- 
Scholes equation. This problem is far more challenging than the one we consider here. 
Indeed, the solution u — u(s, t, a) of equation (12.2|) depends nonlinearly on the coefficient 
o. On the other hand, the problem, which was solved above, is linear. 

The main challenge in solving nonlinear inverse problems is linked with the fact that 
the conventional least squares functionals for them are non convex, thus having many 
local minima and ravines, see, e.g. Isakov (2014) for this observation. Some ideas to 
calculate more accurate values of a ( s,t ) were proposed in Bouchouev and Isakov (1997), 
Bouchouev and Isakov (1999), Bouchouev, Isakov and Valdivia (2002) and Isakov (2014). 
However, these publications have not used our mathematical model. We believe that 
to work with our model, some modifications of currently known globally convergent nu¬ 
merical methods for coefficient inverse problems might be applied, see Klibanov (1997) 
(two papers), Klibanov and Thanh (2014) and Klibanov and Kamburg (2015) for these 
methods. 
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